Prospects for supersymmetric charged Higgs boson discovery 
at the Tevatron and the LHC 
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We investigate the prospects for heavy charged Higgs boson production through the mechanisms 
pp{pp) — > tbH ± +X at the upgraded Fermilab Tevatron and at the upcoming LHC collider at CERN 
respectively. We focus on the MSSM case at high values of tan/3 > mt/mt and include the leading 
SUSY quantum corrections. A detailed study is performed for all important production modes and 
basic background processes for the tibb signature. At the upgraded Tevatron a charged Higgs signal 
is potentially viable in the 220 — 250 GeV range or excluded at 95 %CL up to 300 GeV. At the 
LHC, a H ± of mass up to 800 GeV can be discovered at 5<r or else be excluded up to a mass of 
~ 1.5 TeV. The presence of SUSY quantum effects may highly influence the discovery potential in 
both machines and can typically shift these limits by 200 GeV at the LHC. 
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The full experimental confirmation of the Standard 
Model (SM) is still waiting for the finding of the Higgs 
boson. Last LEP results, suggesting a light neutral Higgs 
of about 115 GeV §, are encouraging, but we will have 
to wait the news from the upgraded Fermilab Tevatron 
or from the upcoming Large Hadron Collider at CERN 
to see this result either confirmed or dismissed. For inter- 
mediate masses above the LEP limit and below 180 GeV 
there is a chance for the Tevatron, but for higher masses 
up to 1 TeV one needs the LHC. However, even if a neu- 
tral Higgs boson is discovered, the principal question will 
stand immutable at the forefront of Elementary Parti- 
cle research: is the minimal SM realized in nature or 
does a model beyond the SM exist with an extended 
Higgs sector? In most of these extensions, the physi- 
cal spectrum contains neutral Higgs particles and some 
of them may mimic the SM one. But in general they 
also involve charged Higgs bosons, and this introduces 
an obvious distinctive feature. For example, in the gen- 
eral two-Higgs-doublet model (2HDM) Q one just adds 
up another doublet of Higgs scalars and then the spec- 
trum of the model contains three neutral Higgs bosons 
h = h°,H°,A° and two charged ones H ± . While the 
detection of a charged Higgs boson would still leave a lot 
of questions unanswered, it would immediately offer (in 
contrast to the detection of a neutral one) indisputable 
evidence of physics beyond the SM. In this Letter we re- 
port on the main results from a fully-fledged study of 
the production in hadron colliders. We restrict our 
analysis to the Higgs sector of the Minimal Supersym- 
metric Standard Model (MSSM), which is beyond doubt 
the most prominent (Type II) 2HDM J|. The lengthy 
details of this work will be presented elsewhere. 

The relevant mechanisms on which we will concentrate 

PPipp) -> tbH* + X (Tevatron) (LHC) , (1) 
are long known to be the leading ones for production 



at high tan/3 ||. They constitute the charged counter- 
part of the process pp — ► tiH+X for Higgs boson produc- 
tion in the SM, recently revisited in j| for the Tevatron, 
and of pp(pp) — * tih + X for neutral MSSM Higgs boson 
production at the Tevatron and the LHC ||. While the 
study of (|I]) has been further addressed in the literature 
to the best of our knowledge all the works on this 
subject- except a first estimation in Q|- do stick to a 
tree-level computation without SUSY quantum effects, in 
spite of the fact that some of them explicitly admit that 
the sort of charged Higgs boson they are dealing with is 
of the MSSM type. Therefore, they are unavoidably af- 
fected by some drawbacks. In the present work we have 
implemented several additional features which improve 
in a substantial manner our knowledge on the real capa- 
bility for the mechanisms (|l]) to produce a charged Higgs 
boson or to put limits on its mass within the MSSM. 
First, we include the leading SUSY radiative corrections 
(both strong and electroweak) along with an analysis of 
the off-shell effects. Second, we perform a beyond-the- 
parton-level simulation of events, which includes the toy- 
detector simulation, jet fragmentation, and initial and fi- 
nal radiation effects. And yet another new aspect of the 
present work is the proper kinematical analysis of the 
gg — > H + tb and gb — > H + i subprocesses at the level of 
differential cross-sections after correctly combining them, 
using the standard recipe for the total cross-section ||, 
and a method based on Ref. |l(| for the differential one. 

A realistic study of charged Higgs boson physics can- 
not be accomplished without including the information 
provided by radiative corrections. These are not only 
potentially large in the computation of the MSSM Higgs 
boson masses themselves but also in the interaction ver- 
tices and self-energies of given processes, particularly in 
the decay of the top quark into charged Higgs and for 
the hadronic decays of the MSSM Higgs bosons 11 
Supersymmetric quantum effects can also be very impor- 
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FIG. 1. Tree-level cross-sections (|lj) versus M H ± at the 
Tevatron Run II (2 TeV) and LHC (14 TeV) for tan/3 = 50. 
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TABLE I. Typical sets of SUSY parameters used in the 
computation of the signal cross-section (0) in Fig. | (all entries 
in GeV). Here \i and M are the higgsino and SU(2) gaugino 
mass parameters, mj is the gluino mass, mj 1 and mj are 
the lightest stop and sbottom masses, and At, At are the top 
and bottom quark trilinear SUSY-breaking couplings. The 
\Sp\ < 0.001 constraint is satisfied. Notation as in Ref. ||llf . 
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FIG. 2. Ksusy = cr SUSY /cr trco factor for the signal <|TJ> as a 
function of tan /3 for the typical sets of parameters in Table^J 
at the Tevatron (for M H ± = 250 GeV) and the LHC (for 
M H ± = 500 GeV). 

tant for the processes (0). Indeed, in Fig. we present 
the tree-level signal cross-sections as a function of Mu± 
for the LHC and Tevatron colliders, whereas in Fig.|| a 
rich variety of SUSY effects is exhibited, as a function of 
tan/3, for the various sets of MSSM parameters indicated 
in Table |. For the signal the leading corrections can be 
described through an effective Lagrangian approach JlJ] 
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gVtb mb tan/3 



V2M W 1 + Am b 
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where rn b is the running bottom mass in the MS. 
The previous formula allows to treat the leading SUSY 
Yukawa coupling effects correctly resummed to all orders 
@. The analytic form of the strong (SUSY-QCD) and 
electroweak (SUSY-EW) corrections Am b in the MSSM 
is given m Ref. |J]. Although Amf, is the only correction 
that contributes at order (a/An) tan/3 (a — as, aw) and 
thus dominates for large tan /3, we have also included off- 
shell SUSY-QCD and SUSY-EW corrections to the iMf* 
vertex and to the fermion propagators. We have made 
extensive use of the CompHEP package for the algebraic 
and numerical calculations p3| . Despite CompHEP is 
only able (in principle) to deal with tree-level calcula- 
tions, with the help of Eq. (||) we have managed to add 
the SUSY corrections to the tbH^ vertex and fermion 



propagators and we have assessed the relevance of the 
off-shell contributions. To this end we have evaluated 
the full set of one-loop SUSY diagrams for the relevant 
tbH + vertex. The same set was considered in detail in 
Ref. jll| for the case where all external particles are on- 
shell. In the present instance, however, at least one of 
the quarks in that vertex is off shell. Therefore, we can 
use the same bunch of diagrams as in the on-shell case 
but we have to account for the off-shell external lines, 
which is a non-trivial task. We have studied this issue 
in detail by expanding the off-shell propagators. First, 
we have modified CompHEP 's Feynman rules to allow 
for the most general off-shell ibH + vertex; then we have 
let CompHEP reckon the squared matrix elements and 
dump the result into REDUCE code. The subsequent nu- 
merical computation of the one-loop integrals has been 
done using the package LoopTools Q. At this point, 
we have inserted expressions for the coefficients of the 
off-shell ibH + vertex that include the one-loop off-shell 
supersymmctric corrections to the vertex itself and to the 
off-shell fermion propagators and fermionic external lines. 
Only half the renormalization of an internal fermion line 
has to be included, the other half being associated to the 
gqq vertex. This procedure has allowed us to estimate 
the relative size of the off-shell effects in the signal cross- 
section, which never exceeds the few per cent level. The 
upshot is that the approximation of neglecting vertex and 
propagator corrections in the cross-section, which may be 
called "improved Born" approximation, is really justified 
in the relevant region of parameter space. 

As for the signal versus background study we have fo- 
cused on the tibb signature corresponding to the H + — > tb 
decay channel, which has the biggest branching ratio at 
high tan/3. Then we have concentrated on the triple- 
b tagging case, which gives the best possibility to mea- 
sure the signal cross-section ([!]). This is a crucial point, 
especially for the Tevatron, where the production rate 
is too small to give any viable signal in the case of 
four-6 tagging. As regards the LHC, a triple-6 tag- 
ging study allows the signal cross-section to be measured 
more precisely, even though the signal/background ra- 
tio can be better for the four-6 tagging case. We no- 
tice that the background processes in Table|H] are m- 
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FIG. 3. Reconstructed tb invariant-mass distribution for 
signal, background and signal plus background events, for 
M H ± = 500 GeV, at the LHC. 

sensitive to the leading type of SUSY effects that con- 
tribute to (0). To perform a realistic signal and back- 
ground event simulation we complied with the following 
procedure. The (tree-level) matrix elements for the sig- 
nal and background processes have been calculated us- 
ing the CompHEP package p3| . The next step was the 
parton-level event simulation, also with the help of Com- 
pHEP. Then we automatically linked the parton-level 
events from CompHEP to the PYTHIA 6.1 Monte Carlo 
generator and the CompHEP-PYTHIA interface [p!s[ . 
Therefore we took into account the effects of the final- 
state radiation, hadronization and string-jet fragmen- 
tation using PYTHIA tools. The following resolutions 
were used for the jet and electron energy smearing: 
AE had /E = 0.8/VE and AE ele /E = 0.2/ VE. In our 
analysis we used the cone algorithm for the jet recon- 
struction with a cone size AR = \J Aip 2 + Ar/ 2 = 0.7. 
The choice of this jet-cone value is related to the crucial 
role of the final-state radiation (FSR), which strongly 
smears the shape of the reconstructed charged Higgs bo- 
son mass. We have checked that the value of 0.7 mini- 
mizes the FSR effects. Now, in order to decide whether 
a charged Higgs cross-section leads to a detectable sig- 
nal, we have to compute the background rate. Since the 
mistagging probability of light quark and gluon jets is ex- 
pected to be < 1% [{l6[ , the only significant backgrounds 
leading to the same ti signature are those shown in Ta- 
ble [II] along with their respective cross-sections. For the 
ttbb and ttqg processes we have applied the jet separa- 
tion cut A]f > 0.5(A fl = yjAO' 2 + A4> 2 ) and the ini- 
tial cut j? T > 10 GeV (Pt > 20 GeV ) at the Tevatron 
(LHC). For the ttj process the initial cut p J T > 10 GeV 
(fy > 20 GeV) was applied at the Tevatron (LHC). To 
obtain a realistic description of the 6-tagging efficiency as 
a function of &-quark transverse momenta, for the Teva- 
tron we use the projected 6-tagging efficiency of the up- 
graded D0 detector || while for the LHC we parame- 
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TABLE II. (a) The main background processes to the sig- 
nal ([[]) at the Tevatron (2nd column) and the LHC (3rd col- 
umn) under the cuts explained in the text. The various contri- 
butions are shown together with that of the subtraction term 
|J; (b) Background from pp — > ttqg when the light quark or 
gluon are misidentified as a 6-jet. 



terize numerical results from the CMS collaboration jl6) . 
Efficiencies for both parameterizations are about 60% at 
the p h T saturation value of ~ 100 GeV. We assume that 
&-jets can be tagged only for pseudorapidity \rjf,\ < 2 by 
both Tevatron and LHC experiments. Furthermore, we 
have optimized the reconstruction procedure, the pr cut 
on the leading 6-jet (j>j, > [M H + /5— 15] GeV) and window 
cut on the ^-invariant mass around the selected values 
of M H + (\m t b — M H +\ < 5y/M H +) to achieve the maxi- 
mal significance of the cross-section signal erg (|l|) versus 
the background, asj\fB. The typical efficiency at the 
Tevatron is 5 — 6% while for LHC it goes down to 1 — 2%. 
These values include the triple ^-tagging, branchings of 
W-bosons decays (leptonic and hadronic decay modes) 
and the efficiency of the kinematical cuts and reconstruc- 
tion of the ttbb signature. As an example, Figure || shows 
the reconstructed tb invariant-mass distribution for sig- 
nal and background events at the LHC. 

In Fig. |] we present the discovery and exclusion limits 
for at the Tevatron and the LHC at high tan/3: It 
shows the signal cross-section as and the cross-sections 
which would lead to the 5a, 3a and 1.98er(95% CL) sig- 
nificance. From the intersections of the last three with 
as we infer the mass ranges which can be discovered or 
excluded. The tree-level case is also shown, and it is seen 
to be too small at the Tevatron to place any limit. But 
when including SUSY effects the situation changes. We 
take the moderate input set B from Table |[ The un- 
known QCD effects at NLO are estimated by including a 
If -factor (typically one expects K ~ 1.5 Q), and we used 
a bottom quark pole mass mj, = 4.6 GeV. For example, 
for M H ± =215 GeV at the Tevatron one would expect 
7 signal and about 6 background reconstructed events 
at L = 25 fb" 1 . At the LHC with M H ± = 500 GeV 
we have 1200 and 3800 signal and background events 
respectively at L = 100 fb _1 . A canonical 5 a discov- 
ery limit around 800 GeV can be obtained at the LHC 
for the MSSM charged Higgs, or else an exclusion limit 
at 95% C.L. up to at least 1.2 TeV (K = 1). For the 
Tevatron we can now place a 95% C.L. exclusion limit 
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SUSY corrections could show up as a smoking gun over 
the tree-level cross-section (after subtracting the conven- 
tional QCD effects) with a contribution of the order or 
larger than the QCD effects themselves. Since both the 
processes (Q) and the SUSY effects are only relevant at 
high tan/3, the sole presence of the signal could be a 
hallmark of the underlying MSSM dynamics. A close 
comparison with the neutral Higgs boson production pro- 
cesses mentioned above should be, if available, very useful 
to confirm or dismiss the MSSM nature of the charged 
and neutral Higgs bosons. 

The work of D.G. and J.G. has been supported by the 
E.U. Marie Curie program, and that of J.S. by CICYT. 
A.B. thanks J. Owens, N. Kidonakis and L. Reina for 
stimulating discussions and acknowledges the support of 
U.S. DOE. 



FIG. 4. Discovery and exclusion limits for the charged 
Higgs boson at the (a) Tevatron and (b) LHC, for tan/3 = 50. 
Shown are the total cross-sections for: i) the signal (lit) at the 
tree-level, ii) the SUSY corrected signal for the input set B 
in Table W including QCD factors K = 1 and K — 1.5, iii) the 
3 a (Tevatron) or 5 a (LHC) discovery limits at the integrated 
luminosities L = J Cdt = 30 fb' 1 and 100 fb' 1 , iv) the 95% 
CL exclusion limit at L = 100 fb" 1 . 



in the mass range 200 — 250 GeV, if K = 1. However, 
if the QCD factor is K ~ 1.5 this would open the ex- 
citing opportunity to observe the charged Higgs already 
at the Tevatron with a 3er significance in the mass range 
220 — 250 GeV!. Moreover, one should notice that b- 
tagging efficiency could be increased up to ~ 70% with 
the use of the 3D vertexing algorithms For triple 

6-tagging, this would augment the signal by at least a 
factor of 2, and so the discovery and the 95% exclusion 
limits would be extended accordingly. Needless to say, 
in case of maximal SUSY enhancement (Cf. set A) the 
exclusion/observation mass range for the charged Higgs 
boson would be further enlarged. 

In summary, we have assessed the possibility to see 
a SUSY charged Higgs at the Tevatron and the LHC 
through the mechanisms (Q). Our study of the quan- 
tum corrections within the MSSM has shown that at 
high tan (3 they are dominated by exceptionally impor- 
tant effects that can be absorbed into an effective tbH^ 
vertex, and therefore in practice they can be treated 
as an "improved Born approximation" . From a de- 
tailed signal versus background study we have shown that 
the prospects for the signal discovery are viable for the 
Tevatron (if M H ± < 250 GeV) and promising for the 
LHC (if M H ± < 1.5 TeV). In deriving these limits we 
have demonstrated that the quantum machinery from 
the MSSM can play a crucial role to increase the sig- 
nal. In large portions of the parameter space the genuine 
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